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ABSTRACT 

We compare [a/Fe], metallicity, and age distributions of globular clusters in elliptical, lenticular, 
and spiral galaxies, which we derive from Lick line index measurements. We find a large number of 
globular clusters in elliptical galaxies that reach significantly higher [a/Fe] values ([a/Fe] >0.5) than 
any clusters in lenticular and spiral galaxies. Most of these extremely a-enriched globular clusters are 
old (t > 8 Gyr) and cover the metallicity range — 1 < [Z/H] <0. A comparison with supernova yield 
models suggests that the progenitor gas clouds of these globular clusters must have been predominantly 
enriched by massive stars (> 20 M©) with little contribution from lower-mass stars. The measured 
[a/Fe] ratios are also consistent with yields of very massive pair-instability supernovae (~ 130 — 
190 Mo). Both scenarios imply that the chemical enrichment of the progenitor gas was completed on 
extremely short timescales of the order of a few Myr. Given the lower [a/Fe] average ratios of the 
diffuse stellar population in early-type galaxies, our results suggest that these extremely a-enhanccd 
globular clusters could be members of the very first generation of star clusters formed, and that their 
formation epochs would predate the formation of the majority of stars in giant early-type galaxies. 
Subject headings: globular clusters: general — galaxies: star clusters — galaxies: formation — galax- 
ies: evolution 



1. INTRODUCTION 

The field of globular cluster system research has con- 
tributed count l ess important insight s to astronomy (e.g. 
Shaplevl 119181 ISearle fe Ziiml 119781 iZinnl 119851 iHarrisI 
1^9lITAsmTLaii^r^rolIT998l and references therein). Be- 
cause they survive a Hubble time, globular cluster sys- 
tems (GCSs) provide important information on the as- 
sembly history of their parent galaxy and the physical 
conditions in the star formation episodes during which 
they formed. With the advent of the Hubble Space Tele- 
scope, extragalactic GCSs in a variety of host galax- 
ies outside the Local Group be came accessible to ac- 
curate photometric studies (e.g. iWhitmore et alJ 119931 
Whitmore fe Schweizerl 119951 iGebhardt fc KisslerT^atid 
19991 IPuzia et alJ 119991 iKundu fc Whitmord l2001albt 
Cote et al.ll2004f) . In the last few years, progress in GCS 
research experienced a boost through the availability 
of efficient multi-object spectrographs, which facilitate 
spectroscopy of large samples of extragalactic globu- 
lar clusters out to large distances (~ 20 Mpc; see e.g. 
Kissler-Patig et al.lll998t iCohen etallll998l IPuzia et all 
20041) . Spectroscopy provides the best constraints on the 
ages and chemical compositions of globular clusters. 

The chemical composition of individual extragalac- 
tic globular clusters provides insight on the enrich- 
ment histories and timescales of their parent massive 
gas clouds. Because of the different progenitor life- 
times of type II and type la supernovae (SNe) (e.g. 
iTornambe fc Matteuccilll986|) an d the different chemica l 
composition of their ejecta (e.g. iNomoto et al.llT997alL 
the progenitor gas clouds will undergo different chemical 
evolution depending on whether they were enriched on 
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short timescales predominantly by type II SNe, or on ex- 
tended timescales by both type II and la SNe. Ages and 
chemical compositions of globular clusters can be used 
to reconstruct the early assembly histories of globular 
cluster systems and their host galaxies, since since self- 
enrichment is excluded for all but the m ost massive glob- 
ular clusters i|Recchi fc Danzigeil 120051) . Together with 
accurate ages and metallicities, the chemical composi- 
tion of globular clusters can be used to reconstruct the 
early assembly histories of globular cluster systems. 

In this paper, we investigate [a/Fe], metallicity, and 
age distributions of individual globular clusters systems 
as a function of environment (in elliptical, lenticular and 
spiral galaxies) based on high-quality spectroscopy data 
from the literature. We identify and discuss extremely 
a-enriched clusters that are (almost) exclusively found 
in ellipticals. 

2. DATA 

The majority of our analysis i s bas ed on the spec- 
troscopic dataset of IPuzia et alJ 1)20041) . who obtained 
high-S/N spectra for globular clusters in seven early- 
type galaxies (NGC 1380, 2434, 3115, 3379, 3585, 
5846, and 7192) with the FORS multi-object spectro- 
graph at ESO's Very Large Telescope (VLT). We aug- 
ment this sample with high-quality index measurements 
from the literature, which a re explained in the last 
section in IPuzia et alJ l|2004fl . We refer the reader to 
this paper for further details. In summary, we add 
to the sample of globular cl usters in ellip tical galax- 
ies data fr om NGC 3610 llStrader et alJ p00¥ and 
NGC 4365 (jLarsen et all 120031 IBroche et a,Ul2005|> . To 
the sample of globular clusters in lenticular galax- 
ies w e add spectra from NGC 3115 llKuntschner et all 
I20021). NGC 4594 l|Larsen et al.ll200l . and NGC 5128 
ijPeng et alJ [2004). The sample of globular clusters 
in spiral galaxies is assembled of clusters in the 
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Milky Wav (IPuzia et alJ 120021) . M31 (iTrager et all 
ITflflg IRfiMlev et 12004 IPuzia et all I2005a|b M81 

(Schroder et alJl2002D. and the Sculp tor spirals NGC 55, 
247, 253, and 300 llOlsen et7al]l200l . 

All index measurements were uniformly per f ormed 
wit h the passband definit i ons of iWorthev et alJ l|1994f) 
and IWorthev fc Ottavianil l)1997|) and are calibrated on 
the Lick standard system using Lick index standards. In 
the following analysis we consider only globular clusters 
which have high-quality measurements in the indices H/3, 
H7A, H(5a, (Fe), Mg&, and Mg 2 since these enter our age, 
metallicity, [a/Fe] fitting routine. To guarantee a robust 
comparison of the three samples, we select data with 
measurement uncertainties AH/3 < 0.5, A H7A < 0-7, 
AH<5 A < 0.7, and A[MgFe]'< 0.4 A (see iPuzia et alJ 
2005b). The final dataset contains 87, 47, and 59 globu- 
lar clusters with high-quality index measurements in el- 
liptical, lenticular, and spiral galaxies, respectively, and 
samples the brightest ~ 10% of the globular cluster lu- 
minosity function. In terms of calibration quality and 
signal-to-noise, it is, to our knowledge, the best-quality 
spectroscopic dataset of extragalactic globular clusters 
currently available. 

3. ANALYSIS 

We use the technique described in IPuzia et al.l (|2005bf) 
to derive ages, metallicities, and [a/Fe] ratios for all 
sample globular clusters. In summary, we compare the 
age-sensitive Balmer line indices H/3, H7A, US a and the 
metallicity-sensitive indices (F e), Mg2 and [MgFel' t o the 
population synthesis models of lThomas et al. (2003) an d 
iThomas et all ()2004|) which provide predictions for stel- 
lar populations with well-defined abundance patterns, 
in particular with varying [a/Fe] ratios. These models 
were previously calibrated on Galactic g lobular clusters 
(jPuzia et alJl2003 : iMaraston et a l. 2003) for which accu- 
rate abundance ratios, including Fe-peak and a-elements, 
were known from high-resolution spectroscopy. Ages, 
metallicities, and [a/Fe] ratios are derived in an iterative 
^-minimization approach of diagnostic grids taking into 
account index measurement errors and systematic uncer- 
tainties of the Lick index system. Since the error distri- 
butions arc similar for all indices that enter the analysis, 
a comparison of all three globular cluster samples with 
the same SSP model predictions guarantees a robust dif- 
ferential analysis. Absolute globular cluster ages have to 
be interpreted cautiously, since the absolute age scale is 
prone to systematic uncertainties in the model calibra- 
tion and an unknown contribution from hot horizontal 
branch stars. However, our analysis allows the distinc- 
tion between old (t>8 Gyr), intermediate-age (2<<<8 
Gyr), and young (t<2 Gyr) globular clusters. The inter- 
nal metallicity and [a/Fe] calibrations are more robust 
and have systematic uncertainties < 0.15 dex. It is im- 
portant to realize that these internal systematics do not 
affect our differential analysis. 

We present [a/Fe] diagnostic plots for globular clus- 
ters in spiral, lenticular, and elliptical galaxies in Fig- 
ure d Figure shows the derived distributions of glob- 
ular cluster ages, metallicities, and [a/Fe] ratios in el- 
liptical, lenticular, and spiral galaxies. The most un- 
biased representation of these data is realized via a 
non-parametric probabilit y density estim ate using an 
Epanechnikov- kernel fsee ISilvermarJ Il986l for details), 
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Fig. 1. — This figure shows [a/Fe] diagnostic plots for globu- 
lar clusters in spiral {upper panel), lenticular {middle panel), and 
elliptical galaxies {bottom panel). Different symbols show different 
globular clusters in individual host galaxies which are described in 
the legend of each panel. Solid lines show iso-[a/Fe] tracks for 0.0, 
0.3, and 0.5 dex for a 13 Gyr old stellar population with metallic- 
ities from [Z/H] = —2.25 to +0.67. To illustrate the systematics 
due to age, we show the same tracks for a 3 Gyr stellar population, 
indicated by dotted lines. Note that the models tracks are virtually 
independent of a ge up to solar metallicities. The model predictions 
were taken from Thomas et al. (2003, 2004). 



which are indicated by solid lines. Dashed lines show 
bootstrapped 90% confidence limits. The median ages for 
the elliptical, lenticular, and spiral samples are 10.2 ±2.1, 
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9.4 ± 2.2, and 10.5 ± 2.0 Gyr 4 . The sampling of all three 
sub-samples is somewhat biased towards the central re- 
gions of galaxies and therefore towards more metal-rich 
globular clusters. In fact, in the elliptical and lenticular 
sample there are only few clusters more metal-poor than 
[Z/H] w —1.5, which correspond to the metal-poor sub- 
popu lation of Galactic halo globular clusters (e.g. iZinnl 
1985). For both the elliptical and lenticular sub-sample 
we find a slightly sub-solar mean metallicity, while the 
spiral sample is much less affected by this bias and probes 
significantly lower globular cluster metallicities down to 
[Z/H] w —2.0. The corresponding median metallicities 
arc [Z/H] = -0.32, -0.23, and -0.58 dex, for the ellipti- 
cal, lenticular, and spiral sample, respectively. 

The perhaps most striking feature in Figure His the 
very broad and possibly bimodal [a/Fe] distribution for 
globular clusters in elliptical galaxies (upper right panel) , 
with a high- [a/Fe] "shoulder" reaching values up to ~ 1.0 
dex. This shoulder is significantly less developed in the 
sample of lenticular galaxies and entirely absent for the 
sample of spiral galaxies. Both spiral and lenticular sam- 
ples have a median value of [a/Fe] is ~ 0.2 dex, twice as 
low as for the sample of ellipticals. We emphasize that, al- 
though NGC 1380 GCS may slightly dominate the high- 
[a/Fe] values in the lenticular sample, no single globular 
cluster system dominates any of the distributions. In- 
stead, globular clusters from all galaxies populate the 
derived distributions rather homogeneously (see Fig. QJ. 

Although there is some evidence for multiple compo- 
nents in the probability density estimate of the [a/Fe] 
distributio n of globular clusters in ellipticals, the KMM 
algorithm (A shman et al.lll994|) rejects the hypothesis 
that the distribution is different from bimodal at a ~ la 
level. The test does not strongly support bimodality, 
but does provide additional evidence that the [a/Fe] 
distribution for globular clusters in ellipticals is signif- 
icantly broader and reaches much higher [a/Fe] values 
than those for globular clusters in lenticular and spiral 
galaxies. 

4. DISCUSSION 

While the averag e properties of the s amples were al- 
ready addressed in iPuzia et al.1 (2005b), we focus here 
on the globular clusters showing extreme [a/Fe] values. 
Figure [3] shows [a/Fe] plotted against age (left panels) 
and metallicity (right panels) for globular clusters in el- 
liptical, lenticular, and spiral galaxies (from top to bot- 
tom). The direct comparison of GCSs shows important 
differences: globular clusters tend to have higher [a/Fe] 
ratios and reach higher metallicities in hosts with earlier 
morphological type. Most of the extremely a-enhanced 
clusters ([a/Fe] > 0.5) are old (t > 8 Gyr), and are 
mainly found in elliptical galaxies. Very few counterparts 
are found in lenticular and spiral galaxies. The metallic- 
ities of these extreme [a/Fe] clusters span a wide range 
from metal-poor [Z/H] rs —2 dex to solar, concentrating 
around values between — 1 < [Z/H] < dex in ellipticals 
and lenticulars, while the few candidates in spirals have 
metallicities [Z/H] < — 1 dex. There is tentative evidence 
for an age- [a/Fe] correlation for globular clusters in el- 
liptical galaxies (correlation coefficient 0.5; see panel a in 
Fig. [SJ, in the sense that younger globular clusters have 

4 The errors give the semi-interquartile ranges. 
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Fig. 2. — Age, metallicity, and [a/Fe] distributions of globular 
clusters in elliptical (upper row), lenticular (middle row), and spiral 
galaxies (bottom row). Solid curves are non-parametric probability 
density estimates with their 90% confidence limits indicated by 
dashed lines. 

on average lower [a/Fe] ratios. For metal- rich globular 
clusters ([Z/H] > —1.0) in elliptical galaxies, we find an 
anti-correlation of [a/Fe] and metallicity, which becomes 
weaker toward later-type hosts. At face value, these re- 
sults suggest that a first generation of globular clusters 
in elliptical galaxies was enriched to much higher [a/Fe] 
values (i.e. faster) than in lenticular and spiral galaxies. 

4.1. Chemical enrichment with different SN types 

We explore in the following whether the [a/Fe] ratios 
of the extremely a-enriched clusters can be reproduced 
by state-of-the-art supernova yield models. For this pur- 
pose we compute [a/Fe] ratios o f the ejecta of type I I 
SNe in the ma ss range 13-70 llNomoto et al.lll997a|) . 
type la SNe l|Nomoto et all Il997b1. hypernovae with 
very large explosion energies > 10 52 erg (jNakamura et alJ 
2001), and pair-instability SNe wit h progenitor masses 
> 160 M Q peger fc WooslevI 12002^ . To stay consistent 
with the definitions of the SSP models, we compute 
[a/Fe] ratios of the SN yields as the mass of the a- 
elements N, O, Mg, Ca, Na, Ne, S, Si, and Ti relative 
to the mass of the Fe-peak elements Cr, Mn, Fe, Co, 
Ni, Cu, and Zn. Th e results are normalized to the solar 
abundances given in lAnders &: Grevessd l|1989ft . We note 
that other combinations of a-elemcnts, in particular the 
addition of carbon, does not affect the results. 

The [a/Fe] ratio of the cumulative ejecta of core- 
collapse SNe with progenitor masses 10 — 50 Mq inte- 
grated over a Salpeter IMF is indicated in Figure 01 as 
a solid horizontal line at ~ 0.36 dex. Such an IMF ex- 
tending to high masses can reproduce the [a/Fe] ratios 
of virtually all globular clusters in spirals and the major- 
ity of globular clusters in lenticular galaxies, but not the 
most extreme [a/Fe] ratios of globular clusters in ellip- 
tical galaxies. A significant contribution of massive stars 
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Fig. 3. — Correlation plots of [a/Fe] ratios with ages (left panels) and metallicities (right panels) for globular clusters in elliptical (upper 
row), lenticular (middle row), and spiral galaxies (bottom row). To illustrate the enrichment by supernovae of different types, mean [a/Fe] 
ratios of the ejecta of type fa SNe (yellow shading), cumulative type II SNe ejecta with progenitor masses in the range 10 — 50Mq (thick 
horizontal line) integrated over a Salpeter IMF, hypernovae with very large explosion energies > 10 52 ergs (cyan s hadin g), and type III 
pair-instability supernovae with progenitor masses 170 — 19OM0 (magenta hatched region) were overplotted (see Sect. UTTl for details). The 
right ordinate indicates progenitor masses of type II SNe that produce the [o/Fe] indicated by the left ordinate. In the upper panels (a 
and b), we plot two enrichment model s for an elliptical galaxy with 10 10 Mq (dashed line) and 10 12 Mq initial baryonic mass (solid line), 
taken from Pipino &z Matteucci 12004). Small crosses indicate average statistical errors. Because of the changing At/t, the statistical age 
uncertainties were calculated for the young and old sub-samples split at 8 Gyr. 



> 20 Mq is necessary to boost the [a/Fe] beyond ~ 0.4 
dex (see Fig.[3J). The ejecta of pair-instability SNe, which 
require very metal-po or progenitor stars, are e xpected to 
be rich in a-elements ijHeeer fc Wooslevl F2002'). Our cal- 
culations indicate that the parent gas clouds of globular 
clusters can be enriched to extreme [a/Fe] ratios by pair- 
instability SNe with masses ~ 130-190 Mq. We find that 
hypernovae with progenitor masses 6 — 16 Mq are ruled 
out as primary contributors to gas with [a/Fe] > 0.4 
since their ejecta are not a-enriched enough. 

We assume that the extreme globular clusters consid- 
ered here have similar masses to Milky Way globular clus- 
ters. E.g., for a typical mass of ~ 10 5 ' 5 Mq, and measured 
metallicities of up to [Z/H] w 0, the total metal content 
of the parent gas cloud amounts to up to ~ 5600 Mq. 
One way of producing this amount of metals is to as- 
sume a progenitor/enriching stellar population composed 
of high-mass (> 20 Mq) stars only (i.e. a progenitor stel- 
lar population with an IMF truncated at the lower end at 
20 Mq), combined with a sufficiently high star- formation 
rate in order to enable the formation of such massive 



stars fsee lGoodwin fc Pagell 2005). In any other scenario, 
the extreme [a/Fe] ratios require extremely short enrich- 
ment timescales, i.e. a very short time between the SN 
explosion of massive stars and the formation of globu- 
lar clusters: of the order of the lifetime of a 20M Q star, 
corresponding to a few Myr. 

Our empirical result that globular clusters with such 
extreme chemical compositions mainly reside in early- 
type galaxies but seem absent in spirals suggests a funda- 
mental difference in enrichment timescales between some 
phase(s) of globular cluster system formation. 

A bimodality in the [a/Fe] distribution, as hinted by 
the top right panel of Fig. [21 would imply a cessation of 
globular cluster formation processes over a period short 
enough to be unresolved in age (on the scale of Figure E3 
~ 10 9 yr), but lasting long enough to allow a significant 
chemical contribution from lower-mass stars (< 20 Mq, 
i.e. few 10 6 yr). Speculating whether such a delay could 
occur in the course of massive galaxy evolution, we note 
that it might be introduced by the reionization epoch 
at redshifts 6 < z < 20, which for a few 10 7 years may 
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have halted, or at least suppressed, the formation of mas- 
sive star clusters and enabled the generation of less a- 
enriched globular clusters which seem to dominate GCSs 
in lenticulars and spirals. 

4.2. Comparison with chemical evolution models 

We compare our observational result for globular clus- 
ters with photochemical evolut i on m odels of elliptical 
galaxies bv lPipino fc Matteuccil l)2004f) . In Figure [3] we 
plot the predictions for [a/Fe] ratios for two galaxies 
with 10 10 and 10 12 M Q initial baryonic mass. The di- 
rect comparison shows that extreme [a/Fe] ratios can 
only be achieved within the first < 0.2 Gyr after the ig- 
nition of star formation processes. Further, the models 
predict similar total metallicities as those of the most 
metal-rich extremely a-enhanced globular cluster in our 
sample. Note that the absolute starburst ignition time is 
arbitrarily set in the models; i.e., one can delay the star- 
burst ignition time by an arbitrary amount which would 
leave the predictions of [a/Fe] as a function of time after 
the burst unchanged. The delay between the bursts of 
10 12 M o and 10 10 M o , on the other hand, is fixed in the 
models, and depends on the total gas mass in a given 
halo. This implies that less-massive starbursts, occur- 
ring later from material with a similar chemical compo- 
sition as for more massive starbursts, can produce glob- 
ular clusters with similarly high [a/Fe] ratios, but with 
younger ages. In a hierarchical galaxy formation picture 
these younger clusters are subsequently merged into more 
massive halos. 

It is important to realize that the very high [a/Fe] 
ratios for old globular clusters in elliptical galaxies are 
not mirrored in the average [a/Fe] of the field stellar 
population as measured from the diffus e light. Using 
the sa me SSP models as in this study, Thom as et al.l 
( 2005) recently found that early-type galaxies have lower 
luminosity-weighted [a/Fe] values (([a/Fe])E ~ 0.25) 
with very few galaxies showing [a/Fe] >0.3 dex, i.e., val- 
ues comparable to those of the low- [a/Fe] peak of glob- 
ular clusters in ellipticals and most clusters in lenticular 
and spiral galaxies. This suggests that the formation 
of massive globular clusters with extremely high [a/Fe] 
ratios predates the formation of the majority of stars 
in elliptical galaxies. However, we emphasize that our 
findings do not exclude the presence of some extremely 
a-enriched field stars in elliptical galaxies. 

Interestingly, this would be consistent with the pop- 



ular scenario in which stars are born in star clusters 
and gradually become part of the diffuse light through 
evaporation by two-body relaxation and tidal stripping 
of th e star clusters (e.g.. lLada fc Ladall2003t [Fall et alJ 
120051) . Indeed, the timescale of such disruption mecha- 
nisms is very long (more than a Hubble time) for massive 
globular clusters of ~ 10 5 ' 5 M Q and shorter for less mas- 
sive ones (M IM oc M -1 / 2 for two-body relaxation; e.g. 
iFall fc Zhand l2001). so that most of the diffuse light in 
ellipticals would originate, in such a scenario, from lower- 
mass clusters, most of which dissolved afte r a few Gyr 
or sooner (see also iGoudfrooii et al.l 12004]) . This actu- 
ally suggests that lower-mass clusters in ellipticals should 
generally not show [a/Fe] > 0.4, which would imply a 
later formation. These predictions should be testable in 
the near future. 

5. SUMMARY 

We present [a/Fe], metallicity, and age distributions of 
globular clusters in elliptical, lenticular, and spiral galax- 
ies, that were derived from Lick line index measurement 
and identify a population of globular clusters with ex- 
treme [a/Fe] ratios. The main results of our study of 
extremely a-enriched globular clusters are: i) they are 
predominantly found in early- type galaxies, ii) most of 
these highly a-enriched globular clusters are old (t > 8 
Gyr) and cover the metallicity range — 1< [Z/H] <0, in) 
a comparison with supernova yield models suggests that 
the progenitor gas clouds of these globular clusters were 
predominantly enriched by massive stars (>20Mq), iv) 
given the lower average [a/Fe] ratios of the diffuse stellar 
population in early- type galaxies, our results suggest that 
the extremely a-enhanced globular clusters are members 
of the very first generation of star clusters formed, and 
that their formation epochs likely predate the formation 
of the majority of field stars in giant early-type galaxies. 
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